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We show that lattice mismatched bilayer graphene can realize robust one-dimensional 
wires. By considering a single domain wall where the masses of the Dirac electrons 
change their sign, we establish a general projection principle. This determines how 
the existence of topological zero-energy domain wall states depends on the direction 
of the domain wall and locations of the massive Dirac cones inside the bulk Brillouin 
zone. We generalize this idea for arbitrary patterns of domain walls, showing that 
the topologically protected states exist only in the presence of an odd number of 
topological domain walls. 
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Two of the most prominent topological phases of matter in condensed matter physics 
occur in the quantum Hall effects^ and topological insulators^. The former exhibits a 
pure topological order- that relies on no unbroken symmetries, while the latter relies on the 
retention of time reversal symmetry. Although a true topological insulator state in graphene 
is difficult to observe because of its weak spin-orbit coupling strength, it is well-established 
that the massive Dirac cones form topological structures - merons - in momentum spaced 
Therefore topologically protected states can still emerge at the domain walls and edges of 
gapped graphene and bilayer graphene 6 -"— . However, these states are yet to be realised in 
the lab. 

The purpose of the current work is two-fold. Firstly, we propose an ideal material for 
experimentally realizing one- dimensional wires as topological domain walls: mismatched 
substrate or bilayer graphene. The canonical example of lattice-mismatched graphene is a 
twisted bilayer. This material is 'trivial' in the context of our results. Two 'topological' 
examples are bilayer graphene where one layer has a single line defect^, or where one layer 
is uni-axially strained. The most topical example of mismatched substrates at present, is 
graphene atop hexagonal boron nitride 13 . 

The parallel purpose of this work is to understand the general properties of multiple 
domain walls in samples where the phase-relaxation length is large in comparison to the 
length scale of the domain wall patterns. We find that the prescription of the conditions 
under which robust zero-energy domain wall states will exist, contains two ingredients. First, 
we establish a general projection principle, which determines how the existence of topological 
zero-energy domain wall states depends on the direction of the domain wall and locations 
of the massive Dirac cones inside the bulk Brillouin zone. Then, we generalize this idea for 
arbitrary patterns of domain walls and show that the topologically protected zero-energy 
states exist only in the presence of an odd number of topological domain walls, resulting in 
a novel type of even-odd effect. 

It is well known that at the interface between two massive Dirac systems with opposite 
signs of mass, which we shall call a domain wall, a zero-energy localized state exists. This 
observation forms the basis for the existence of surface states in topological insulators, and 
has also been pointed out in the context of gapped graphene 6 . Creating an oscillating mass 
term is not trivial. In this work we induce a mass on an effective single layer Dirac system 
by spatially varying the inter-layer coupling in a bilayer system. The connection between 
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FIG. 1. The two necessary ingredients to realize robust one- dimensional wires, (a) Spatially 
varying mass: in mismatched bilayers, interlayer coupling maps onto an effective Dirac mass, (b) 
The generalized projection principle: in the particular case domain walls in graphene, the projected 
cones overlap and so may become gapped in the armchair case, whereas for zig-zag walls they remain 
distinct. 

on-site energy and Dirac mass has been established elsewhere^, and the principle is outlined 
in Fig. Ufa). We assume that the states with different spins do not couple. The effective 
theory of a single spin species system with N massive Dirac cones is 



where, for example in graphene there are two cones described by du2) = (±vp(k x — 

1(2) 1(2) 

K x ),VF{ky — Ky ),m), where Vp = 3t/2 is the Fermi velocity, and m is a constant 
mass. The low energy theory for spatially varying mass is obtained with the replacements 
tti y 77i(r), and k x{y) -)■ -id x(jj) . 

Firstly, we introduce a projection principle which defines the topological nature of the 
domain wall states for a bulk material with any number of Dirac cones. For simplicity, 
we consider first a single domain wall and assume that the low-energy spectrum of the 
bulk material is described with two Dirac cones as shown in Fig. []Jb). Because of the 
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translational invariance of the system along the domain wall, the component of k in that 
direction remains a good quantum number and the states in the 2D Brillouin zone of the 
bulk material become effectively projected onto a ID Brillouin zone, potentially causing 
mixing of the states with the same value of k along the domain wall direction. Therefore, 
we can expect two possible regimes to arise from the projections depending on the domain 
wall directions, as schematically illustrated in Fig. [1] (b). The first regime, which we shall 
call the topological regime, is that where the low-energy parts of the two cones are projected 
onto different regions in the ID Brillouin zone. The second regime, which we call the trivial 
regime, is that where the projection causes the two valleys to overlap at low energies. It turns 
out that this distinction between topological and trivial regimes is necessary to determine 
the existence of robust zero-energy eigenstates localised on the domain wall. Similarly, for 
a multiple Dirac cone bulk system, the domain wall states caused by a given Dirac cone are 
topologically protected if the projected Dirac cone retains the distinction from any of the 
other projected Dirac cones. 

Algebraically applying this principle to Hamiltonian (JTJ, we consider a domain wall at 
y = 0, and look for solutions where E = 0, obtaining ip^ x (x,y) = exp[ik x x + iK y y — 
Vp 1 Jq dy'm(y > )) (1 — 1) per cone. If the two cones are at positions K and K', as in graphene, 
then the two projected solutions are separated by K x — K' x . In the trivial regime however, the 
two low energy valleys overlap, and admixing of the zero energy solutions occurs, analogous 
to the case of an armchair ribbon^. In this case then, the zero energy solutions split into 
two particle-hole symmetric finite energy solutions. 

Projecting these solutions onto a ID domain wall parallel to d, we obtain the effective 
theory for an N cone bulk system 

N 

H(k) = J2s i (k-K i -d) (2) 

i 

where Sj = ±1, depending on the corresponding dj in the Hamiltonian ([1]). This is a 
generalization of the field theory proposed in reference^, and is simply a Luttinger liquid. 
Importantly, in the absence of scattering of magnitude |(K* — K J ') • d\ for all j, we have a 
single distinct zero-energy surface solution arising from cone i. 

Having established the projection principle for the presence of a number of surface states, 
we now introduce the notion of an even-odd effect for domain wall states: A coupled 2n — 1 
topological domain-wall system has at least one single zero-energy mode localized on the 
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FIG. 2. Three mismatched bilayer graphene systems, (a) shows a twisted bilayer that hosts 
topologically trivial domain walls. The walls that meet in the centre of the figure are shown by 
a dotted line. In (b) the upper layer has been uni-axially strained, thus allowing a set of equally 
spaced vertical topological domain walls. Sample (c) contains a single line defect, where a line of 
blue atoms has been removed from the top layer, resulting in a step-like change of mass along the 
topological domain wall. 

domain walls. This definition can be seen to be true from a simple bandstructure argument. 
For 2n — 1 domain walls, the mass changes sign from, say, positive to negative n times, and 
negative to positive n — 1 times. This leads to n solutions at the K point with positive 
slope, and n — 1 solutions at the K point with negative slope. When these solutions are 
admixed, we must do so pair-wise. Each pair becomes gapped, and we are left with a single 
zero energy solution at the K point with positive slope. Similarly we are left with a single 
zero energy solution at the K' point with negative slope. 

Further, we can quantify this argument as follows. Consider a system of Nw domain walls 
described by Eq. (j2J), where the distance between neighbouring walls is d, and the overlap 
integral of these nearest neighbors is to- The Hamiltonian around the point k x = K x is 
given by 



N w 



H = ^v F {k x - K)cl xd c kxJ + ^t D {c\ x j+l c kx j + c.c. 



(3) 



where tp is the overlap integral of states on neighboring walls. In this case, the solution to 
the inter-domain wall coupling in the Hamiltonian takes the tri-diagonal form 




FIG. 3. Explicit bandstructure calculation demonstrating the even-odd effect in a typical zig-zag 
graphene nanoribbon under an oscillating staggered sublattice potential. N\y is the number of 
domain walls in the sample, and we clearly see that when N\y/2 € Z only, the system is gapless. 



H' - diag JV (e) -> 



-e t D 

t*D —£ 
t D -e 



\ 



V 



(4) 



7 



Searching for e = solutions, we find that only odd dimensional matrices support such 
solutions. This confirms the even-odd effect. We now turn to the discussion of lattice- 
mismatched bilayer graphene. 

It has been shown recently that twisted bilayer graphene displaying a Moire pattern has 
an oscillating effective mass 14 . This is a direct manifestation of the effect introduced in 
Fig. HJa). In order to develop some intuition for this, consider two single, uncoupled layers 
of graphene. If we introduce only the dominant interlayer tunnelling term 7, then in the 
normal Bernal-stacked regime, 7 directly couples the A sublattice of the upper species with 
the B sublattice of the lower, as the upper A site sits directly above the B site. However, 
if the A-sublattices constitute interlayer nearest-neighbours, then 7 couples them instead. 
If we construct a single-layer effective theory by integrating out one of the two layers, then 
this interlayer coupling becomes an effective mass. The effective Hamiltonian is then none 
other than Eq. ([1]). We should emphasize that lattice- mismatching of single layer graphene 
with a substrate is in principle completely equivalent. 



The pattern of domain walls in a twisted bilayer exhibiting a Moire-like pattern is shown 
in Fig. |2^a). We have numerically confirmed that in the most promising situation of small 
rotation angles, the domain walls are oriented along the armchair direction. Two systems 
with 'topological' walls are shown in Fig. Efb) and (c). The former is a set of parallel domain 
walls, obtained by uni-axial straining of the upper layer. This results in a continuously 
varying mismatch between the two lattices, and thus again a continuously varying, oscillating 
mass term. Fig. [2]Jc) shows a bilayer system where the upper layer has a single line defect. 
A line defect is obtained by removing a single line of A or B atoms, thus translating the 
entire lattice a small constant amount. Rather than a continuous variation of mass, this 
system contains a sharp step-like change in sign of the mass. For systems (a) and (b), the 
number of domain walls in the sample is tunable by either varying the angle of rotation or 
the amount of uni-axial straining, respectively. 

Finally, we confirm numerically the existence of a single zero-energy state in systems with 
an odd number of topological domain walls. In Fig. [3] we show the dispersions between two 
projected K points for typical zig-zag nanoribbons with Nw = 1,2,3 and 4 domain walls. We 
find that in the presence of an odd number of topological domain walls only a single domain- 
wall mode is gapless and all the other modes have an energy gap which is determined by the 
tunneling energy. This energy gap can be significant if the domain walls are close to each 
other. We note that we have also confirmed the triviality of armchair domain wall systems, 
thus demonstrating the generalized projection principle. 

In conclusion, we have shown by general topological arguments, that topologically triv- 
ial insulators with multiple massive Dirac cones, where a local change in mass induces a 
domain wall, support robust topological states. By a generalized 'projection principle', we 
showed that the topological regime is that where the projection of the Dirac cones onto the 
domain wall does not cause the two cones to overlap. We also defined an even-odd effect for 
these systems, which demands that a system with an odd number of coupled domain walls 
necessarily supports at least a single zero energy mode. Mismatched bilayer graphene was 
then shown to be an ideal system to experimentally realize such states. 
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